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ABSTRACT 

We present a three epoch survey for transient and variables in the extended Chandra 
Deep Field South at 5.5 GHz with the Australia Telescope Compact Array. A region 
covering ~0.3 deg^ was observed on timescales of 2.5 months and 2.5 years and typical 
sensitivities 12.1 — 17.1 pJy beam“^ (Icr) were achieved. This survey represents the 
deepest search for transient and variable radio sources at 5.5 GHz. In total 124 sources 
were detected above the 5.5a level. One highly variable radio source was found with 
AS > 50% implying a surface density of ^3 deg“^. A further three radio sources 
were found with lower levels of variability equating to a surface density of ^^13 deg“^ 
above a detection threshold of 82.3 p,Jy. All of the variable sources have inverted radio 
spectra (between 1.4 and 5.5 GHz) and are associated with active galactic nuclei. We 
conclude that these variables are young gigahertz peaked-spectrum sources with active 
and self-absorbed radio jets. We explore the variability completeness of this sample 
and conclude that the fairly low levels of variability would only be detectable in 3—25% 
of all sources within the field. No radio transients were detected in this survey and we 
place an upper limit on the surface density of transient events < 7.5 deg“^ above a 
detection threshold of 68.8 ^Jy. 

Key words: instrumentation: interferometers, radio continuum: general, techniques: 
image processing, catalogues 


1 INTRODUCTION 


There is a vast amount of parameter space available to 
search for transient and variable radio sources. The physical 
mechanisms driving dynamic radio emission remain poorly 
constrained as a function of frequency and flux density. This 
is especially true at higher frequencies (e.g. 5 GHz and 
higher) where the field of view (and survey speed) of ra¬ 
dio interferometers decreases. At a range of survey cadences 
and at frequencies around 1.4 GHz, the dyn amic sky has 
been well characterised above ^ one mjy fe.g. iBo wer et abl 
2010l: Croft et alll201Cll : [Bower et al.ll20lj] : ICroft et al. 20 111 : 
Croft et al. I l2013li . At 5 GHz and below one mJy however, 
this is not the case. 

At 1.4 GHz one of the benchmark surveys in explorin g 
the sub-mjy variable radio sky is that of ICarilli et al.l EoOSli . 
In this survey, five archival repeat VLA observations of the 
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Lockman Hole were examined for variable radio sources on 
timescales of 17 months and 19 days. Of the sources detected 
between 50 — 100 /iJy, 2% were found to be variable. The 
surface density of variable sources above lOO^Jy was found 
to be ~ 18 d eg~^. On the shorter timescales of one day to 
three months. I Moolev et al.l (l2013l i examined archival obser- 
vations of the exten ded Chandra Deep Field South (eCDFS; 
iLehmer et akll^OOSl l. At 1.4 GHz a small fraction (1%) of 
sources were found to vary significantly above 40 /iJy. 

At ~5 GHz, the deepe st variability survey to date is 
that of iBecker et al.l 1I2OIOI '). which probed a flux density 
greater than one mJy, an order of magnitude higher in flux 
than t hat of the Carilli et al. (2003) survey. The lBecker et all 
(l2010li survey showed that the surface density of variable 
sources towards the Galactic Plane above one mJy is ~1.6 
deg“^ on timescales of days to years. In comparison Gre¬ 
gory & Taylor (1986) showed that a much stricter limit of 
~ 10~® deg“^ is found above 70 mJy on timescales of days 
and years. To date there has been no blind survey at 5 GHz 
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that has studied the abundance of long duration variable 
radio sources below 1 mjy. 

A few blind radio surveys have searched for both tran- 
sien ts and variab l es at 5 GHz on long timescales. For exam¬ 
ple, |B^^7^eF^ (l2007l ~l (also see lFrail et al.l[2oH ~) searched 
archival data (at 5 and 8 GHz) of a single calibration held 
observed over 22 years with the Very Large Array (VLA). 
Upper limits wer e placed on t he su rface density of long du¬ 
ration transients. lOfek etliD (1201 il l searched for both tran¬ 
sients and variables using the VLA at 5 GHz. One possible 
radio transient was reported and 2% of the source sample 
(> 0.5 mJy) were found to vary on a timescale of two years. 

Here we present a new set of observations aimed at 
probing the time-domain properties of the radio sky at 
5 GHz below 1 mJy. The eCDFS is a ~0.3 deg^ region 
favoured for its low optical and HI extinction. A plethora 
of multi-wavelength data from both space and ground based 
facilities are available for this region (e.g. Grogin et al. 2011, 
Xue et al. 2011, Koekemoer et al. 2011, Skelton et al. 2014). 
At 1.4 GHz, a 3.6 deg^ region encompassing the eCDFS 
has been observed by the Australian Large Array Survey 
fATLAS: iNorris et al.lf200^ team using the Australia Tele¬ 
scope Gompact Array (ATGA). At both 1.4 GHz and 5 GHz, 
observati ons of the eCDFS have also been m ade using the 
VLA ('see lKellermann et al.l200^ . In late 2009. [Huvnh et ^ 
(I 2 OI 2 II obtained observations of the eCDFS at 5.5 and 9 GHz 
(~ 144 hours). 

For the survey presented in this paper we obtained 
a further two epochs of observations (~100 hours) of the 
eCDF S at 5.5 GHz (designed to match those of iHuvnh et al.l 
I2OI2I I. The motivation for this survey was to search a be¬ 
spoke region of parameter space (i.e. 5 GHz and <1 mJy) 
for signatures of dynamic radio emission on long timescales: 
months to years. It was also to survey a region which had ex¬ 
tensive multi-wavelength coverage, thus to negate the need 
for potentially lengthly follow-up observations. The multi¬ 
wavelength data allow for accurate source identification, 
which has impeded the physical interpretation of discoveries 
in previous surveys. 

The outline of this paper is as follows. In sections 2 
and 3 we describe the observations and data reduction proce¬ 
dure, respectively. In section 4 we describe the transient and 
variable search methodology. Sections 5 and 6 describe the 
source counts and characteristics of the variable and tran¬ 
sient sources detected, respectively. In section 7 we discuss 
the physical interpretation of the results and we compare 
our work with that reported in the literature. 


2 OBSERVATIONS 

For the analysis presented in this paper we used archival 
data of the eGDFS held (PI: M. Huynh, project code C2028, 
see iHuvnh et al.ll2012l l and new observations, both of which 
were obtained with the ATGA. We will refer to the archival 
data as epoch 1 henceforth. These archival observations were 
conducted between 2009 August 12 and 2010 January 16 
and we include all data, regardless of the date collected, 
into this epoch. In total 115 hours hours of data were col¬ 
lected and the observations were typically organised into 12 
hour blocks. Approximately 80% of the observations were 
obtained in 2010 January. The hnal day of observations for 


epoch 1 were obtained on 2010 January 16 and we take this 
as the reference time for examining the variability. 

We obtained two further epochs of data for this sur¬ 
vey (project code G2670). In Table [T] the dates, integration 
times and array configurations of the observations are sum¬ 
marised. The observations for epoch 2 were performed from 
2012 May 31 and 2012 June 03 . The observations for epoch 3 
were obtained between 2012 August 14 to 2012 August 19. 
In total we sampled timescales between epoch 1 and epoch 2 
of 869 days (or 2 years, 4 months and 18 days) and between 
epoch-2 and epoch-3 of 77 days. This was calculated between 
the final dates of the observations for each epoch. 

We ad o pted the same observing strategy as 
iHuvnh et al.l (I 2 OI 2 II to reproduce, as best as possible, 
the observations of epoch 1. We used 42 pointings in a 
hexagonal layout to mosaic a field of view ~ 30' x 30'. The 
42 pointings were spaced such that at 5.5 GHz the map 
was sampled to uniform sensitivity. We observed each of 
the 42 pointings consecutively (interleaved with calibration 
observations) over the course o f the observing block (as 
was done by iHuvnh et al.ll2012l L For epochs 2 and 3 the 
observations were organised into five ~10 hour blocks and 
executed intermittently over the dates specified above. We 
obtained approximately 48 hours of observations on epochs 
2 and 3, which resulted in a ~ decrease in sensitivity 
for these epochs, when compared with epoch 1. 

All of the images have been restored with a restoring 
beam of 5.0" x 2.0" with beam parallactic angle (BPA) 1.4 
degrees. We used the full available continuum bandwidth 
of 2048 MHz centred at 5.5 GHz using t he Gompact Array 
Broadband Backend JWilson et al.ll201ll L The observations 
were centred at RA 3^32"'2® and Dec —27°28'34" (J2000). 


3 DATA REDUCTION 

There have been some improvements to the Multichan¬ 
nel Image Reconstruction, Image Analysis and Display 
fMIRIAD: ISault et al.|[l9^ ~l software package and the wide¬ 
band imaging steps used to reduce the data since the anal¬ 
ysis of Huynh et al. (2012). We follow the data reduction 
methodology described by Huynh et al. (in prep) but sum¬ 
marise it here briefly. Epoch 1 from Huynh et al. (2012) was 
re-processed to be consistent with the latter epochs. Ampli¬ 
tude, phase and bandpass calibration solutions were derived 
from intermittent observations of PKS B1934-638 (ampli¬ 
tude and bandpass calibration) and PKS J0348-2749 (phase 
calibration). 

For each epoch of observations, each of the 42 pointings 
was deconvolved, GLEANed, and self-calibrated separately. 
The final image for each epoch was produced from a linear 
mosaic of the separate pointings. In Figure [T] we show an 
image of epoch 2. In epochs 1, 2 and 3 we achieved noise 
levels of 12.1, 17.1 and 16.2 ^Jy beam“^ respectively, close 
to the centre of the image (see also Table [1}. 

Figure [ 2 ] shows the noise profile of the three epochs as a 
function of sky area. For each epoch we generated an r.m.s 
map, where each map pixel represented the Icr r.m.s sen¬ 
sitivity in that region. For each epoch we summed up the 
cumulative area (deg^) and plotted this as a function of sen¬ 
sitivity. Figure [2] shows that because the epochs are made 
of linear mosaics the noise profile is relatively flat out to 
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Table 1. A summary of observations used for this survey. Epoch 1 was obtained by Huynh et al. (2012). In total we have probed 
timescales of 869 days (between epochs 1 and 2) and 77 days (between epochs 2 to 3). The N^et column gives the number of blind source 
detections in each image. The RMS column gives the typical noise calculated close to the centre of the image. 


Epoch 

Date 

Total integration time 
(hours) 

Array 

RMS 

(pjy beam“^) 

^Det 

1 

2009 Jun 08 

- 2010 Jan 16 

115 

6A, 6D 

12.1 

115 

2 

2012 May 31 

- 2012 Jun 03 

48 

6D 

17.1 

82 

3 

2012 Aug 14 

- 2012 Aug 19 

48 

6A 

16.2 

91 
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Figure 1. An image of epoch 2. A total of <^0.3 deg^ was imaged using 42 separate pointings. In the central region of the map an r.m.s 
sensitivity of 17.1 /rjy beam“^ was achieved. A total of 48 hours of data were used to make this image and 82 sources were detected. 
See Table [T] for details of other epochs. All images have been restored with a restoring beam of 5.0” X 2.0” with beam parallactic angle 
(BPA) 1.4 degrees. 


~0.2 deg^. Table [2] shows the mean sensitivity over five cu¬ 
mulative sky area bins for epoch 1. These values are used in 
subsequent analysis to define the upper limit on flux density 
and area that this survey probed. 


4 TRANSIENT AND VARIABLE SEARCH 
TECHNIQUES 


An automated pipeline developed as part of the Aus- 
tralian Square Kil ometre Array Pathfinder (ASKAP; 
|johnstoiiet_gl SOO^JVariable and Slow Transients (VAST) 
project ( Murphy et al ] l2013ll was used to search the obser¬ 
vations f or transient and v ariable radio sources. We refer the 
reader to lBell et al.l (l2014[ l for a full description of the tran- 
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Figure 2. Sensitivity <Jcut (/^Jy bearn”^) versus cumulative area 
A (deg^) as a function of epoch. The cumulative area gives the 
amount of sky we have observed with sensitivity cr < acut- For 
example, in epoch 1 approximately 0.4 deg^ is observed with sen¬ 
sitivity better than 100/iJy. 




( 2 ) 


The values of xL expected to follow the theoretical 
distribution Xt (for n — 1 degrees of freedom), assuming 
no variability is present. Using the cumulative distribution 
function (CDF) we therefore calculate the probability P that 
the xL value could be produced by chance. A variable source 
was defined as having P < 0.001. This sets a false detection 
rate at 0.1 per 124 sources searched (total source counts 
discussed below). 

Secondly we calculated the de-biased modulation index 
for all source light-curves (see Akritas fc Bershadvl Il996l : 


iBarvainis et al.l l2005l : lSadle£^^jl 


index (see equation 5 in iBell et al. 


200^. The modulation 


2014 ) does not take into 


account the errors on the flux measurements. A bias there¬ 
fore exists whereby sources with large errors (or low signal to 
noise) have larger modulation indices. The de-biased modu¬ 
lation index corrects for this effect. The de-biased modula¬ 
tion index is defined as 


rrid = 


1 ELi(s.-5)2-Er=i<T? 


(3) 


Table 2. The mean sensitivity of epoch 1 calculated over five 
discrete sky area bins. 


Cumulative Area 
(deg2) 

Mean Sensitivity 
(/ijy bearn”^) 

0.1 

11.7 

0.2 

12.5 

0.3 

15.5 

0.4 

26.0 

0.5 

63.0 


sient and variable search methodology. To summarise briefly, 
we located sour ces within the imag es using the AEGEAN 
source finder fsee lHancock et ^l2012l ~l. We used a source ex¬ 
traction level of 5.5(7 to search for sources. A forced measure¬ 
ment algorithm was used to produce complete light-curves 
for all sources detected in one image greater than 5.5 (j. For 
all the light-curves (constructed from either blind detec¬ 
tions or blind detections plus forced measurements) we then 
tested for significant variability using the following statistics. 

Firstly, we tested for variability by calculating the x^ 
probability that the source remained constant. The mean of 
the data points were used as a model to test if the light-curve 
was invariant. For each source we calculated xL where 


2 _ ■sp (Si - Sf 


( 1 ) 


Si is the ith flux density measurement with variance af.nis 
the total number of epochs. The weighted mean flux density, 
S, is defined as 


where S is the mean of the flux density measurements. 
When md is imaginary the magnitude of the variability is 
considered uncertain , due to the size of the errors Ui (see 
IBarvainis et al.ll2005l l. If rud is imaginary we evaluate the 
modulus of the function inside the square root. We then 
multiply the final expression by minus one i.e. 


md = -lx 


1 



n 


(4) 


This allows us to compute and plot all value s f or md - 

For comparison with ICarilli et al.l (l2003l i and 
iMooIev et al.l l|2ni3l i we also calculate the fractional 
variability, defined as: 


AS 


Smax Sn 


s 


(5) 


where Smax and Smin are the maximum and minimum 
flux density measurements for a given source. We classify 
a highly variable source as having AS > 50%. 

In this dataset we classify a variable as a source that 
was detected in all epochs and showed significant variability 
(based on the x^ statistic). We classify a transient source as 
one that was detected in less than three epochs, with large 
signal to noise ratio. The distinction is fairly trivial but the 
physical processes driving the different classes may possibly 
be different e.g. explosive and one off vs. intrinsic and persis¬ 
tent. In conjunction with the variability statistics described 
above we inspected all sources that were not detected in all 
three epochs. 
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5 SOURCE COUNTS 

In total 124 Gaussian source components were detected in all 
three images. Table [T] gives the number of blind detections 
as a function of epoch. Depending on the image and noise 
level different Gaussian components were blindly detected 
in extended objects. A total of 115 sources were blindly de¬ 
tected in the deepest image (epoch 1). A further nine Gaus¬ 
sian components were detected in epochs 2 and 3 (and not 
epoch 1) due to the complex morphology of extended ob¬ 
jects. 

In the analysis of iHuvnh et al.l ll2012h 142 source com¬ 
ponents are reported in epoch 1. The discrepancy between 
the source counts in epoch 1 (115 vs. 142) can be attributed 
to different source finding algorithms and source extraction 
levels. The VAST pipeline represents each source as a com¬ 
bination of Gaussian components. Unresolved or partially 
resolved sources are characterised as single elliptical Gaus- 
sians, whilst resolved sources are characterised by a number 
of such components. Over all three epochs, for each blindly 
detected Gaussian fit, a complete light-curve was produced 
(using forced measurements if the source was undetected in 
one or more epochs). These light-curves were then tested for 
significant variability (as described in section 4). 


6 RESULTS 

6.1 Variable sources 

In Figure |3] we plot the variability statistics as a function of 
signal-to-noise (SNR) ratio for all of the sources in the sam¬ 
ple. The large amount of sources with negative de-biased 
modulation indices is a result of having two epochs that are 
less sensitive. The source finder is very efficient at detecting 
objects in epoch 1, which are subsequently measured in the 
less sensitive epochs (2, 3) with slightly larger error bars. To 
check that this was expected we wrote a simulation. Using a 
Gaussian number generator we simulated the light curves of 
124 sources which had error bars in proportion to those from 
the three epochs. For this simulation the de-biased modula¬ 
tion index shifted to be predominantly negative for sources 
with low SNR. We then repeated the calculation, this time 
keeping the error bars (for the three data points) of equal 
size. This yielded an even distribution of modulation indices 
around zero. 

In total seven sources were identified which met our 
criteria of being variable. After visual inspection one of 
the sources was found to be part of an extended multi- 
component object and we discarded the variability because 
the Gaussian fits were confused, extended and unreliable be¬ 
tween images. Two further sources were discarded because 
the morphologies were extended and the de-biased modu¬ 
lation indices were low. Four unresolved sources were con¬ 
sidered to have significant astrophysical variability and we 
will discuss the individual properties of these below, also 
see Figure |4l The properties of these variable sources are 
summarised in Table (3] 

6.1.1 ATCDFS ,7033209-274249 

This source showed the largest de-biased modulation index 
md = 24.4% (see Figure |4] upper left panel) and fractional 
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Figure 3. A plot of the de-biased modulation index versus 
all sources in the sample. The size of the marker is proportional to 
the mean signal-to-noise ratio of the source detection. The dashed 
line denotes the variability detection threshold. Seven sources sit 
above the variability detection threshold. Four of these sources 
are considered genuine variables and are used in the subsequent 
analysis. 


variability A5 = 57.9% and it was the faintest of the four 
variab les. This source was not identified as being variable 
in the iMoolev et al.l (l2013l l sample. It has a spectroscopic 
redshift of j* = 0.733 (Szokoly et al. 2004) and is classified 
as a type-one AGN from its X-ray properties (Afonso et al. 
2006). The HST AGS F180LP image from CANDELS (Gro- 
gin et al. 2011, Koekemoer et al. 2011, Skelton et al. 2014) 
shows a spheroidal optical host with a nearby companion to 
the south west (see Figure [S] top left). This companion has 
a photometric redshift of 0.73 (Ferreras et al. 2005). The 
variable source and its companion appear to be members of 
a cluster at 2 = 0.734 (Gilli et al 2003, Szokoly et al. 2004, 
Silverman et al. 2010). 

This source has an inverted radio spectrum with 
= -1-0.52 where S oc (Huynh et al. in prep). This 
spectral index was calculat ed between archival measure¬ 
ments bv iMiller et al.l ll2013ll at 1.4 GHz and measurements 
in a deep image of all three epochs at 5.5 GHz (Huynh et al. 
in prep). The spectral indices between 1.4 GHz and 5.5 GHz 
for the individual epochs are ai = -1-0.65, 0:2 = -1-0.29 and 
a = -1-0.29. The intra-band spectral index is relatively flat 
until 5.8 GHz, above which point it becomes negative (see 
Figure [ 5 ]). The intra-band spectral index is measured in a 
deep image of all three epochs (Huynh et al. in prep). 


6.1.2 ATCDFS J033211-273726 

This source had the highest xic value (lowest P-value) of all 
the variables and a de-biased modulation index = 9.2% 
(see Figure |4] upper right panel). A significant dip in flux 
density is observed in epoch-2 with respect to epochs 1 and 
3. The source had previously been identihed bv iHuvnh et al.l 
ll2012ll to be variable with r espect to archival measurements 
by Kellermann et aP ll2008h and we confirm that continued 
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ATCDFS J033209-274249 



Time 



Time 


Figure 4. Light-curves of the variable sources listed in Table |3] F 
indicated in the plot. The top left light-curve has a y-scale range 
range of ±25%. 


ATCDFS J033211-273726 



Time 



Time 


each of the sources, the xf^ de-biased modulation index are 
’ ±50% of the mean flux. The remaining light-curves have a y-scale 


Table 3. Properties of the variable sources. The SNR column denotes the mean signal-to-noise ratio of the flux measurements. The 
“l 4GW2 column gives the spectral index (where S oc measured between the 1.4 GHz flux density (taken from [filler et ahll2013ll 

and 5.5 GHz averaged data from all three epochs (Huynh et al. in prep). See Huynh et al. (2012) for positional errors on individual 
source positions. 


Name 

RA (J2000) 

Dec (J2000) 

SNR 


p-value 

rrid (%) 

AS (%) 

^5.5GHz 

^lAGHz 

Host 

ATGDFS J033209-274249 

3h32m9.7s 

-27d42m48.4s 

29.3 

18.3 

1.0 X 10"'^ 

24.4 

57.9 

±0.53 

AGN 

ATGDFS J033211-273726 

3h32mll.7s 

-27d37m26.3s 

629.2 

368.5 

< 1 X 10-15 

9.2 

20.9 

±0.83 

QSO±AGN 

ATGDFS J033208-274734 

3h32m8.6s 

-27d47m34.6s 

217.3 

44.5 

2.2 X 10-1° 

7.0 

17.9 

±0.41 

QSO±AGN 

ATGDFS J033206-273235 

3h32m6.1s 

-27d32m35.8s 

428.4 

14.3 

7.8 X 10-1 

1.9 

4.9 

±0.15 

AGN 


radio variability via this survey. This source was also iden - 
tified as being variable at 1.4 GHz bv iMoolev et al.l ll2013l) . 
This source has an R-band (Vega) magnitude of 18.626 and 
is cla ssified as having a QSO-like optical SED (jWolf et al.l 
l2008h . An optical spectrum of the source shows broad emis¬ 
sion lines indicative o f quasar activity and p laces the source 
at redshift 2 = 1.575 llSilverman et al.ll201(]|) . 

The radio spectral index is af, 4 GHz = ±0.83 (Huynh et 
al. in prep), so it is an inverted spectrum source. The intra¬ 
band spectral index is = ±0.75, see Figure (5] The 

broad optical emission lines and inverted radio spectrum 
suggests a line of sight relatively face-on to the torus of the 
AGN with a view to the accretion disk. This is an X-ray 


luminous source with L(0. 2 - 2 keV) = 2 x 10'*^ erg s“^, 
classifying it as an AGN (IXue et al.l l201l| ). Figure [6] (top 
right panel) shows the co mpact optical co unterpart in the 
HST ACS F180LP image (iRix et al.ll2004^ . confirming the 
QSO nature. 

6.1.3 A TCDFS ,7033208-274 734 

The light-curve of this source (Figure |4] bottom left panel) 
shows an increase in flux over all three epochs with nid = 
7.8%. This source w as not identified as being variable in the 
iMoolev et al.l (l2013l ) sample. It has been classified as a QSO 
using the ultraviolet excess method (Groom et al. 2001) and 
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Figure 5. Intra-band spectral energy distributions for the four 
variable sources. These spectra are measured in a deep image of 
all epochs (Huynh et al. in prep) across the 2 GHz observing 
bandwidth which was centred at 5.5 GHz. 


Table 4. Upper limits on the surface density of transient sources 
on timescales of 2.5 years and 2.5 months, as a function of detec¬ 
tion threshold. The detection thresholds are defined as the mean 
sensitivities given in Table [2] multiplied by the source extraction 
level of 5.5 ((t). 


Cumulative Area 
(deg^) 

Detection threshold 
(/ijy beam“^) 

Surface density p 
(deg-2) 

0.1 

64.4 

<15.0 

0.2 

68.8 

<7.5 

0.3 

82.3 

<5.0 

0.4 

143.0 

<3.8 

0.5 

346.5 

<3.0 


6.2 Transient sources 

No transient radio sonrces were fonnd above 5.5 ct. As dis¬ 
cussed in section 4, by using onr forced measurement algo¬ 
rithm and inspecting light-curves with less than three blind 
detections we would have been capable of detecting any such 
transient source above this limit within this dataset. 


it also has a QSO-like optical SED (Wolf et al. 2008). The 
optical spectrum shows emission lines indicative of an AGN 
at z = 0.544 (Groom et al. 2001). This source is classified 
as a type-one QSO from its X-ray luminosity and hardness 
ratio (Szokoly et al. 2004). 

The source has an inverted radio spectrum ai\ 4 aHz = 
-1-0.41 (Huynh et al. in prep). It has an intra-band spectral 
index of = +0.69, see Figure[5l This source is point¬ 

like in HST ACS F180LP image from CANDELS (Grogin 
et al. 2011, Koekemoer et al. 2011, Skelton et al. 2014) (see 
Figure [6l bottom l eft). This sour c e was also identified as an 
optical variable bv lTrevese et al.l (l2008l ) and potentially the 
radio and op t ical v ariability are intrinsically linked e.g. see 
iBreedt et al.l (l2010l ). 


6.1.4 ATCDFS ,1033206-273235 

This source showed a very low level of variability rrid = 1.9% 
(see Figure 0] bottom right panel) and was the brightest of 
the variable sources and all other sour ces in the field . It was 
not identified as being variable in the lMoolev et all (l2013l ) 
sample, but the NVSS 1.4 GHz flux density differs from the 
later epoch eCDFS observations by 3 mjy (almost 40%). 
Regardless of the low de-biased modulation index the vari¬ 
ability is considered significant and also the radio spectral 
index is intriguing, so we include it in our variables sample. 

It has a reliable spectroscopic redshift of 2 = 0.961 (Sil¬ 
verman el. al. 2010), with the spectrum only showing narrow 
emission lines. It is formally unresolved in our 5.5 GHz radio 
imaging but marginally resolved in higher resolution VLA 

1.4 GHz imaging (Miller et al. 2013). Its X-ray luminosity 
L(0.2 - 2 keV) of 3 x 10^® erg s“^ (Vattakunnel et al. 2012) 
places it in the regime of X-ray luminous AGN (e.g. Bauer et 
al. 2004). This source has an inverted radio spectrum with 
<^i'. 4 GHz ~ +0.15 and intra-band spectral index of a^^an^z 
— +0.40, see Figure (5] 


7 DISCUSSION 
7.1 Transient sources 

No radio transients were detected in this survey. We can 
therefore place upper limits on the surface density of events. 
Using Poisson statistics fsee lGehrelslll98^ with zero detec¬ 
tions at the 95% confidence level we can set an upper limits 
on the surface density (p) of events per deg^ as: 


ln(0.05) 

(n-1) X H’ 


( 6 ) 


where n is the number of epochs and Q. is the sky area 
surveyed (deg^). We use the sky area and corresponding 
mean sensitivities defined in Tabled] to calculate the upper 
limits on the surface density given in Table 01 Within a 
region of 0.2 deg^, where the sensitivity profiles of the 
images are fairly flat (see Figure 01), we find an upper limit 
p < 7.5 deg“^ above 68.8 p4y on timescales of 2.5 months 
and 2.5 years at 5.5 GHz. 


These calculated surface densities are dependent on the 
transient evolution timescale. The cadence of these observa¬ 
tions are not explicitly folded into equation l6l but they are 
implied by n i.e the number of epochs that are separated 
by a given cadence, or in this case, two different timescales 
(2.5 years and 2.5 months). For incohe rent sources the typ¬ 
ical evolution tim escale is < 100 days ilMetzger et al.ll2015l : 
IPietka et al.ll20f5h . which is well matched by our 2.5 month 
timescale (or 77 days). The quoted surface densities of tran¬ 
sients therefore may be even lower for long duration incoher¬ 
ent sources (e.g. 2.5 years), as the long duration sampling is 
not well matched to the evolution timescale (for this exam¬ 
ple). 
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-3>’32'”09.5* -3''32"'07.0* 

Right Ascension (J2000) Right Ascension (J2000) 

Figure 6. Optical HST images of the variable sources ATCDFS J033209-274249 (top left), ATCDFS J033211-273726 (top right) and 
ATCDFS J033208-274734 (bottom left). These panels are HST ACS F180LP images from CANDEL S (Gr ogin et al. 2011, Koekemoer et 
al. 2011, Skelton et al. 2014). The bottom right panel shows a MUSYC i band image dCawiser et al.ir2006l ') of ATCDFS J033206-273235. 
The radio contours from this work start at five times the image noise and increase by multiples of two. 


7.2 Comparison with previous transient surveys 

In comparison to previous work, iBower et al.l (l2007l ) report 
no radio transients on the timescale of one year, over 17 indi¬ 
vidual epochs, with a typical flux density threshold of 90 /rJy 
(at the half-power point). This equates to an upper limit of 
<2.1 deg“^ at the 95% conhdence level (using equation [6] 
assuming a held of view of 0.09 deg^ and 17 epochs). This 
upper limit is shown in Figure [7] and is more constraining 
than our upper limits at comparable detection levels. 

iBell et ^ (1201 ll ) searched for transients in 24 years of 
archival Very Large Array at 1.4, 4.8 and 8.4 GHz. An upper 
limit of p < 0.032 deg“^ was placed on the surface density 
of tra nsients > 8 mjy with typical timescales of 4.3 to 45.3 
days. iBell et ahl ll201ll) were sensitive to shorter timescale 
transient events at a higher detection threshold over a range 


of frequencies. It is therefore difficult to c ompare the surve y 
presented in this paper with the results of iBell et al.l ll201ll) . 

lOfek et’ahl (l201ll) report the possible detection of a 
transient at 5 GHz with a peak flux density of 2.4 mJy. The 
transient was detected in only one epoch at the start of the 
survey, which sampled timescales of days to years. Because 
the transient was detected i n the first epoch o f the survey 
the duration was uncertain. lOfek etld] (1201 ll ) report that 
a sky surface density of transients 0.039 deg“^ is expected 
above 1.8 mJy. 

Scaling this upper limit to a detection threshold of 
0.1 mJy (assuming a euclidean population with p oc 
where Sv is the detection threshold) we would expect to find 
p < 3.0 deg“^, or 0.89 sources over our held of view. This is 
just below the upper limits on the surface density of tran¬ 
sients placed through this work. It is therefore unlikely that 
we would have been able to detect a transient of this type. 
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T.3 Variable sources 

We find four variable radio sources on timescales of months 
to years. With a total of 124 sources this equates to 3% of 
the sample that were variable. All of the variable sources are 
contained within a region covering 0.3 deg^. Due to the high 
signal to noise ratio of the source detections and relatively 
flat noise profile of the images, all four variables would have 
been detectable within any part of the 0.3 deg^ region. This 
gives a surface density of 13.3 variable sources per deg . 

All four of the variable sources had positive spec¬ 
tral indices between the average 5.5 GHz measurements 
and 1.4 GHz measurements taken in 2007 by Miller et al. 
(2013). ATCDFS J033211-273726 had a spectral index of 
a = -1-0.83. For this source over year long timescales we 
measure a fraction variability of 20.9% (see table [3]). Ad¬ 
justing the average 5.5 GHz flux density accordingly we can 
place errors on the spectral index of a = 0.831 q ^3 on this 
timescale. Therefore if the current level of 5.5 GHz variabil¬ 
ity is typical, it is difficult to reconcile this source as having 
been flat spectrum, or even negative in the recent past (as¬ 
suming the 1.4 GHz flux has remained relatively constant). 
The intra-band spectral indices confirm the inverted spectra 
for three of the variable sources (see Figure [S]). 

The inverted radio spectra combined with the opti¬ 
cal and X-ray identifications suggest that these are com¬ 
pact sources, that are possibly young and have active radio 
jets. Synchrotron self absorption can cause inverted radio 
spectra during outburs ts, or shocks, in relativistic jets (see 
iMarscher fc GeaJill985l l. Typically the outburst is followed 
by an optically thin phase resulting in a negative spectral 
index (at the reference frequency of 5.5 GHz). The optically 
thin phase is a function of frequency and it will evolve from 
high to low frequencies as a function of time. This model 
seems most applicable to ATCDFS J033209-274249 i.e. the 
intra-band spectral index is fairly flat, and we could inter¬ 
pret the light-curve and inverted archival spectral index to 
suggest a recent outburst. The spectrum of this source also 
evolves from -1-0.65 to +0.29. This supports the shock-in-jet 
model because the higher frequency end of the intra-band 
spectrum is negative and the averaged spectral index per 
epoch is evolving to become negative . 

Alternatively these are young radio sources with giga¬ 
hertz peaked-spectrum (GPS). The variability could there¬ 
fore be related to relatively recent jet activity. Variability 
has been reported before in GPS sources by, for example, 
iBolton et al.l (l2006l l. In this study they note that the ma¬ 
jority of sources that had their spectra peak below 5 GHz 
were non-variable, whereas almost all the sources which had 
their peak above 5 GHz were variable. The variable sources 
ATCDFS J033211-273726, ATCDFS J033208-274734 and 
ATCDFS J033206-273235 have peaks in their spectra above 
5 GHz and may well be explained by GPS variability. 

The abundance of GPS sources is difficult to quantify 
for a sample at a given flux level and frequency. Typically 
the turnover in a gigahertz peaked source spectrum ca n 
vary fr om 500 MHz up to 10 GHz or more llO’Deairi998ll . 
lO’Deal (119981 ) report that the expected abundance of GPS 
sources withi n a sam ple greater than 1 Jy is ~10% at 5 GHz. 
iHuvnh et al.l (120121 ') report that of the 108 eCDFS sources 
(non multi-component) that had reliable 1.4 GHz cross¬ 
matches, 10 (or 9% of the sample) had spectra af;! > +0.5. 


This compares with 3% of the sources that we find to be 
variable and have inverted spectral indices. Clearly not all 
of the steep spectrum sources are variable within our sam¬ 
ple, however, this assumes that we were capable of detecting 
the variability in the weaker GPS sources (see discussion in 
section 7.5). Given the expected abundance of GPS sources 
in this sample, these variables do constitute a fairly rare 
region of the spectral population. 

We crosschecked the seven vari a bles r eported in the full 
variability analysis of iMoolev et al.l (l2013l) with our sample. 
ATCDFS J033211-273726 (discussed above) was the only 
source found to be vari able in both surveys . The difference 
between this survey and iMoolev et akl (l2013l 'l was the typical 
timescale o f the observat i ons a nd the frequency (1.4 GHz 
vs. 5 GHz'). iMoolev et al.l ll2013h sampled typical timescales 
of one day to three months whereas this survey samples 
~2.5 months and ~2.5 years. The timescales for relativistic 
jets to evolve at these frequencies will be longer than the 
scintillation timescale and hence will be better selected by 
our observin g strategy (i.e . time scale and frequency) when 
com pared to Moole y et al.l (|2013[ ). 

[Moolev et_alj __l[20lj) used the Miller et al. (2013) 
1.4 GHz and iHuvnh et al.l (l2012l ~l 5.5 GHz catalogs to cal¬ 
culate the spectral indices. They report that four variables 
had negative spectral index, and three were positive (with a 
range —0.4 to +0.89). Flat spectrum sources are often com¬ 
pact and can undergo refractive scintillation, therefore the 
difference in spectra l properties between the variables in our 
sample, and those in iMoolev et al.l (l2013l) . may be influenced 
by the presence of scintillation. 


7.4 Comparison with previous variability surveys 

In Figure [7] we compare the surface density of variables 
foun d in this su r vey w ith those reported in the litera¬ 
ture. ICarilli et al.l (l2003l ) report a surface density of variable 
sources above lOO/rJy to be <18 deg~^ at 1.4 GHz on 17 
month timescales (out to a radius of 7.8'). We find a com¬ 
parable surface density ~ 13.3 deg“^, if all the variables 
sources within our sample are considered and we ignore the 
change in frequency. This calculation is performed by sim¬ 
ply dividing the number of variables found (four) by the area 
(0.3 deg"). 

If we assume that our resu lts are d r awn fr om a Poisson 
process, using the method of iGehreld lll986l ) we can cal¬ 
culate upper and lower limits on the surface density. For 
our results, assuming four detections with a survey area of 
0.3 deg " over three epoch s we find 2.3 < p < 15.3 deg". 
For the lCarilli et al.l (l2003l ) results, if we assume that eight 
variables were detected over three epochs with an area of 
7r(32')" = 0.89 deg" we find 2.2 < p < 8.1 deg". Ou r resul ts 
are therefore in broad agreement with ICarilli et all (l2003l l. 

We find only one source (ATCDFS J033209-274249.0) 
that is considered highly variable with AS/S > 50% (see Ta¬ 
ble [3]). This gives a surface density of highly variable sources 
3.3 deg“" (or 1% of our sample). Using Poisson statistics this 
equates to 0.1 < p < 8.0 deg". 

The predicted number of variable sources is of course a 
func tion of flux d e nsity, cadence, field of view and frequency. 
The lCarilli et al] (l2003f ) survey had some similar properties 
to our survey i.e. (i) rms noise measurements (in the range 
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Figure 7. Surface density (deg~^) of transient and variable sources versus flux density (Jy) for a selection of surveys. This plot is 
focused primarily on surveys at 5 GHz but also includes relevant surveys at other frequencies. Superscript V and T indicate variable and 
transient surface densities respectively. CDFS^ gives the surface density of 13 deg“^ at a flux density > 82.3 /rjy for all the variables 
detected in this paper. CDF^ gives the upper limits on surface density of transients calculated through this work. The dashed grey 
line shows a fit normalised to the CDFS'^ surface density where p oc Si, This assumes that the variable source population is drawn 
from a homogenous euclidean population. The additional surveys are labelled accordingly. Gre86^ - Gregory & Taylor (1986); Car03^ 
- Garilli et al. (2003); FraOS^ - Frail et al. (2003); BowOT^ - Bower et al. (2007); BeclO^ - Becker et al. (2010); Ofell^ - Ofek et al. 
(2011); Belll^ - Bell et al. (2011); Mool3^ - Mooley et al. (2013). 


15—18 /iJy beam ^); (ii) cadence (timescales of both months 
and years); and (iii) held of view (~ 0.5 deg^). 

At low Galactic latitudes, lOfek et’ahl (l201ll ) reported 
that > 30% of unresolved sources brighter than 1.8 mjy were 
significantly varia ble on the timescales of days (at 5 GHz). 
lOfek et"^ ll 201 ll) conclude that the likely cause for the short 
timescale variability is refractive scintilla tion. On the longer 
timescales of two years iCTek et al.l | 20 l 3 ) report that a much 
lower fraction (^2%) of the sam ple is variable (>0.5 mJy). 
The findings of lOfek et al.l (1201 ll ) on long timescales are con¬ 
sistent with the results presented in this paper. 

_ On the timescales of days to 15 years, iBecker et al.l 

ll 2010 l) examined the variability of sources towards the 
Galactic plane in the flux density range 1—100 mJy. A sur¬ 
face density of 1.6 deg~^ varia bles was reported ( 8 % of 
the sample). iBecker et al.l (120101) argue that 80% of the ob- 
j ects i n their sample were Galactic objects. The IBecker et all 
(I2OIOI) sample could therefore differ greatly from the objects 
in our sample, otherwise this is a signihcant increase. 

iGregorv fc Tavl^ (Il986l) find a lower surface density of 
variables towards the Galac tic plane. Above ~ 7 5 mJy on the 
timescales of days to years IGregorv fc Tavl^ (119861) report 
~ 2% of the sample was variable at the > 50% level (surface 
densit y < 10“^). This is four times less than IBecker et ^ 
1I2OIOI) although at a higher flux density cut off. 


T .5 Variability completeness 

To assess our ability to detect variability within this dataset 
as a function of source signal-to-noise ratio and variability 
amount and type, we calculate a variability completeness 
metric. We take the light-curves of the four variables sources 
and re-calculate a hypothetical value after scaling the 
source flux density measurements by a factor 7 where 0 < 
7 < 1. By reducing the flux density of the measurements (of 
the variable source light-curves) the variability becomes less 
significant when compared with the errors, as a function of 
7 . The Xc variable is defined as: 


2 _ ilSi - Sc)^ 


( 7 ) 


The modified weighted mean flux density. Sc, is defined as 



( 8 ) 


where 7 is the scaling term. 

Figure | 8 ] shows the values as a function of 7 
for the four variable sources. From Figure |5] it can be 
seen that ATCDFS 033211-273726, which had the high¬ 
est value of xfc = 368, can be reduced in flux by 81% 
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Figure 8. The hypothetical Xc values calculated for the four vari¬ 
able sources as a function of 7 . The 7 factor scales the flux density 
measurements of the variable source light-curves (see equation[7]|. 
The dashed black line shows the variability detection threshold 
(x^ >13.9) used in this survey. 

(7 = 0.19) before the variability is non-detectable. Con¬ 
versely ATCDFS 033206-27273235, which had the lowest 
value of xfc = 14.3 and smallest de-biased modulation index 
rrid = 1.9%, could only be reduced by 2% (7 = 0.98). 

We can translate the point (yum) at which the variabil¬ 
ity becomes non-detectable into a mean signal-to-noise ratio 
SNRvar, for each of the four variable sources. We then use 
SNR^ar to calculate how many other sources within the 
field have a SNR greater than SNR^ar- The results are 
summarised in Table [S] 

Interestingly the fairly lows levels of variability re¬ 
ported in this survey would only be detectable in a small 
number of sources. At best, the variability observed in 
ATCDFS J033209-274249 could be detected in 31 other 
sources (or 25% of the sample) within the field. The variabil¬ 
ity seen in ATCDFS J033206-273235 could only be detected 
in four other sources (or 3% of the sample). This is assum¬ 
ing that there are 124 sources within the field and some 
of these were multi-component or extended objects, which 
are unlikely to be variable. We can of course rule out 100% 
variations in the majority of sources. Our statistics are in¬ 
fluenced by the small number of epochs we have obtained 
(three). Acquiring more epochs should make the statistic 
more robust to detecting significant variability. 


8 CONCLUSION 

We have surveyed the extended Chandra Deep Field South 
for transient and variable radio sources. This survey searches 
a deep (< 1 mjy) part of time-domain parameter space 
on long timescales at 5.5 GHz. We found four sources that 
showed significant variability, for one of which the variability 
was at a very high level. No radio transients were detected 
and we place upper limits on the occurrence of such events. 


Through this survey we have explored the spectral and 
multi-wavelength properties of variables within the eCDFS. 
All of the variable sources had positive spectral indices and 
were associated with QSO and/or AGN. We conclude that 
the physical interpretation for one of the variables is best 
explained by episodic jet activity common to AGN. We con¬ 
clude that three of the variables are young GPS sources 
whereby the variability and spectra are indicative of fairly 
recent activity within the radio jet. 

The radio spectra of these sources have implications for 
future variability surveys. Inverted spectra sources are quite 
rare within a typical radio survey and by definition, these 
sources will be more difficult to detect at lower frequencies 
(e.g. 1.4 GHz). Assuming that the variables discovered in 
this paper are drawn from a much larger population (over 
the whole sky), long duration variability surveys may be 
more efficient targeting higher frequencies to search for such 
objects. 

Papers describing variability surveys often quote the 
number of variables detected as a fraction/percentage of 
the total number of sources. This is misleading because it 
does not factor in the ability to detect different types and 
amounts of variability within those objects. Furthermore, 
the detection threshold is typically used as a proxy for the 
threshold at which variability could be detected. This is 
true, however a source at the detection threshold of a sur¬ 
vey would have to have large and significant variability to 
be detected. Typically (at 5 GHz) we are dealing with fairly 
low levels of variability (from AGN activity) and this is only 
detectable in the brightest sample of sources. 

Within this dataset we have achieved sensitivities be¬ 
tween 12.1 to 17.1 /iJy beam“^ (l<^)i utilising ~250 hours of 
telescope time. We can clearly rule out extremely large am¬ 
plitude variations in the bulk of sources. However, even in 
optimistic scenarios we are probably only capable of prob¬ 
ing low levels of AGN type variability in sources >~ 1 mJy. 
This is approximately two orders of magnitude higher than 
the sensitivity of this survey. It is also worth noting that the 
variables we have detected (which have inverted spectrum) 
con stitute a fairly rar e type of AGN (~10% of the sam¬ 
ple; [Huynh|eU^[20l3). The requirement to be both greater 
than ~1 mJy (for variability to be detected) and belong to 
a class of object in the minority (i.e. GPS), means that we 
have detected variability in the majority of inverted spec¬ 
trum sources that we were capable of detecting it within. 
We could therefore from a statistical standpoint conclude 
that variability is prevalent in these types of objects. 

Detections of transients still remain rare. Large area, 
sensitive, high cadence observations are now required to 
explore regions of par ameter space not co vered by previ¬ 
ous surveys (e.g. see iMetzger et al.l [201^. A number of 
next generation wide-field GHz radio interferometers that 
meet these requirements will soon survey the sky for tran¬ 
sient and variable radio source s. These include: the Ka¬ 
roo Array Telescope (MeerKAT: iBooth et al.ll2009l) and the 
Australian Square K ilometre Array Pathfinder (ASKAP; 
Ijohnston et al.ll2 008l) a nd potentially the Square Kilometre 
Array (SKA: iGarilli &: Rawlingjl2004h . A comparative sur¬ 
vey to this one over 1000 deg^ with, for example, MeerKAT 
could detect may thousands of these variable sources. If 
an unknown transient population exists, leveraging the in- 
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Table 5. Table summarising our ability to detect the variability found in the four variable sources within all sources in the field. SNR 
is the mean signal-to-noise ratio of the variable detections, 'yum is the limit on the flux scaling parameter below which the variability 
becomes non-detectable. SNR^ar is the mean signal-to-noise ratio at 'yum- The Nsources > SNR^ar column shows how many other 
sources within the field have a mean signal-to-noise ratio greater than SNR^ar- 


Source 

SNR 

'yiim 

SNRvar 

^sources ^ SNRvar 

ATCDFS J033211-273726 

629.2 

0.19 

119.5 

12 

ATCDFS J033208-274734 

217.3 

0.55 

119.5 

12 

ATCDFS J033206-273235 

428.4 

0.87 

372.7 

4 

ATCDFS J033209-274249 

29.3 

0.98 

28.7 

31 


creased survey speeds of next generation instruments will be 
the optimum path to discovery. 
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